Obtaining the Neutrino Mixing Matrix with the Tetrahedral Group 



A. Zee 

Kavli Institute for Theoretical Physics 
University of California 
Santa Barbara, CA 93106 
USA 

zee@kitp.ucsb.edu 
(to appear in Physics Letters B) 

We discuss various "minimalist" schemes to derive the neutrino mixing matrix using the tetra- 
hedral group Ai. 

I. Neutrino mixing matrix 

The neutrino mixing matrix V relates the neutrino current eigenstates (denoted hy a = e, /i, r, and coupled by 
the W bosons to the corresponding charged leptons) to the neutrino mass eigenstates (denoted by t-i, i = 1, 2, 3 and 
endowed with definite masses m^) according to 
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Thanks to heroic experimental efforts, the neutrino mixing angles have now been determined [1] to be given by 
sin^ 9i2 - 0.31, sin^ 02^ ~ 0.50, and sin^ 6*31 ^ ( 
(with C23 = cos 023, S23 = sin 023, and so forth) 



sin^ 9i2 ~ 0.31, sin^ 023 ~ 0.50, and sin^ ^31 ~ 0.01, with the mixing angles defined by the standard parametrization 
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(This parametrization may differ slightly from others in that we take detV23 = det V12 = —1.) The error bars are 
such that 031 is consistent with 0, in which case the CP violating phase e"^ does not enter. 

We could suppose either that the entries in V represent a bunch of meaningless numbers possibly varying from village 
to village in the multiverse landscape as advocated by some theorists of great sophistication or that they point to 
some deeper structure or symmetry as some theorists with a more traditional faith in the power of theoretical physics 
might dare to hope for. It is natural to imagine that there is a family symmetry [2] linking the three lepton families. 

Starting with the standard model we assign (all fermionic fields are left handed) the lepton doublets tpa ^ ^ 

the lepton singlets {a — 1,2,3), and the required Higgs fields to various representations of a family group [3] Gp- 
Indeed, if we guess that S12 = l/\/3, S23 = l/\/2, and S31 = 0, we obtain the attractive mixing matrix 
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first proposed by Harrison, Perkins and Scott [4]. Later, X. G. He and I independently arrived at the same Ansatz 
[5] Also, this mixing matrix (but curiously, with the first and second column interchanged) was first suggested by 
Wolfenstein more than 20 years ago [6] based on some considerations involving the permutation group S3. It has 
subsequently been studied extensively by Harrison, Perkins and Scott [7] , and by Xing [8]. Attempts to derive this 
mixing matrix have been discussed by Low and Volkas [9] [10]. A parametrization of the experimental data in terms 
of deviation from V is given in [11]. Following Wolfenstein and defining Vx = {>^^^ + ^'t)/\/2 and Uy = (t'^ — i't-)/^/^, 
we see that (5) says that the mass eigenstates are given by 



1 



Z^l = - Y ^t'e + "yl^^' ^2 = + Y t'S = yy (6) 

The basis {i^i, vi) is rotated from {i^e, v^^ through arcsin(l/\/3) ^ 35°. 

In this paper we will take the neutrinos to be Majorana [12] as seems likely, so that we have in the Lagrangian 
the mass term £ = —VaMo^jiCvp + h.c. where C denotes the charge conjugation matrix. Thus, the neutrino mass 
matrix M is symmetric. Also, for the sake of simplicity we will assume CP conservation so that M is real. With this 
simplification, the orthogonal transformation V^MV produces a diagonal matrix with diagonal elements mi , m2 and 
ms- We are free to multiply V on the right by some diagonal matrix whose diagonal entries are equal to ±1. This 
merely multiplies each of the columns in V by an arbitrary sign. Various possible phases have been discussed in detail 
in the literature [13,14]. 

At present, we have no understanding of the neutrino masses just as we have no understanding of the charged 
lepton and quark masses. The well-known solar and atmospheric neutrino experiments have determined respectively 
that Am| = m2 - mf ~ 8 x 10~^ eV'^ and Aml^.^ = ml-ml ^ ±2.4 x 10~^ eV"^. The sign of Aml^.^ is currently 
unknown, while Amg has to be positive in order for the Mikheycv-Smirnov-Wolfenstein resonance to occur inside 
the sun. We could have either the so-called normal hierarchy in which [mal > |to2| ~ |wi| or the inverted hierarchy 

|"^3| < |»Tl2| ^ 

II. Family symmetry and the tetrahedral group 

For some years, E. Ma [15] has advocated choosing the discrete group A4, namely the symmetry group of the 
tetrahedron, as G_f. With various collaborators he has written a number of interesting papers [16] [17] [18] [19] using 
A4 to study the lepton sector. 

For the convenience of the reader and to set the notation, we give a concise review of the relevant group theory. 
Evidently, A4 is a subgroup of SO{3) (which was often used in the early literature on family symmetry but which 
has proved to be too restrictive.) Since the tetrahedron lives in 3-dimensional space, A4 has a natural 3-dimensional 
represcintation denoted by 3 suggestive of the 3 families observed in nature. The tetrahedron has 4 vertices and 
thus A4 is also formed by the even permutations of 4 objects so that A4 has 4!/2 = 12 elements which could 
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be represented as elements of SO{3). Besides the identity 7 = | 1 ] ,we have the 3 rotations through 180° 

1 

, r2 = I 10 I , and r"3 = I —10. Then we have the cyclic permutation c = 






which together with ricri,r2cr2 and rscrs, form an equivalence class with 4 members. Finally we have 
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the anticyclic permutation a = | 1 | , which together with riari,r2ar2 and r^ars, form another equivalence 

1 

class with 4 members. Thus, the 12 elements belong to 4 equivalence classes with membership 1, 3, 4, and 4, which 
tells us that there arc 4 irreducible representations with dimension dj such that = 12 which has the unique 

solution di = d2 = ds = 1 and ^4 = 3. The natural 3-dimensional representation 3 has just been displayed explicitly. 

The multiplication of representations is easy to work out by using the following trick. Start with the familiar 
multiplication within SO{3) :3x3=l-|-3-|-5. Given two vectors x and y of 50(3), the 3 is of course given by the 

cross product xx y while the 5 is composed of the symmetric combinations 2:22/3 -|- 2:31/2, x^yi + xiy^, Xiy2 + X2yi, 
together with the 2 diagonal traceless combinations 2xiyi — X2y2 — x^y^ and X2y2 — x^y^. Upon restriction of 5*0(3) 
to A4 the 5 evidently decompose into 5 ^ 3 -|- 1' -|- 1" with the 3 given by the 3 symmetric combinations just displayed. 
The 1' and 1" could be taken respectively as linear combinations of the 2 traceless combinations just given: 

1' r^u' = xiyi + ujX2y2 + oj'^xays (7) 

and 

1" ~ u" = xiyi + uJ^X2y2 + ujX3y3 (8) 
with w = e'^'^/^ the cube root of unity so that 
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+ = 0. (9) 

It is perhaps worth emphasizing the obvious, that while 1' and ]_" furnish 1-dimensional representations of A/^ they 
are not invariant under A^. For example, under the cyclic permutation c, u' uju' and u" uP"ul' . Evidently V x 
1" = 1, 1' X 1' = l",and 1" X l" = l', and also (l')* = l". 

Thus, under A/^ we have 3 x 3 = 1+ 1_' + 1" + 3 + 3. It is perhaps also worth remarking that the two 3's on the right 
hand side may be taken as (2:22/3, xaj/i, 2:12/2) and (2:32/2,2:12/3,2:22/1). The existence of 3 inequivalent 1-dimensional 
representations also suggests the relevance of A^ to the family problem. I can't resist mentioning here the possibly 
physically irrelevant fact that [20] alone among all the alternating groups A„'s the group A^ is not simple. 

III. A minimalist framework 

Given these attractive features of ^4, there has been, perhaps not surprisingly, a number of recent attempts [15] 
[21] [22] to derive V using A4. In our opinion, they all appear to involve a rather elaborate framework, for example 

supersymmctry, higher dimensional spacetimc, and so on. Within this recent literature Ma [27] has produced a 
particularly interesting and relatively economical scheme in which the neutrino mixing matrix depends on a parameter 
such that when that parameter takes on "reasonable" values the matrix V as given in (5) is recovered approximately. 

The guiding philosophy of this paper is that we would like to have as minimal a theoretical framework as possible. 

Within a minimalist framework, charged lepton masses are generated by the dimension 4 operator 

O4 = V^l'^'^l^ (10) 

Here denotes generically the standard Higgs doublet, of which we may have more than one. According to a general 
low energy effective field theory analysis [23] [24] [25] neutrino masses are generated by the dimension 5 operator 

O5 = (Cr2^)C(rT2^) (11) 

in the Lagrangian. Here ^ and ^' denote various Higgs doublets that may or may not be the same as the (/?'s. We will 
suppress the charge conjugation matrix C and the Pauli matrix T2 in what follows. It is important to emphasize that 
the analysis leading up to (11) is completely general and depends only on SU{2) x U{1), and not on which dynamical 
model you believe in, be it the seesaw mechanism or some other mechanism (such as the model in [26]). 

We suppose that the family symmetry remains unbroken down to the scale of SU{2) x U{1) breaking, so that the 
operators O4 and O5 have to be singlets under Gp. As is completely standard, when (p, ^, and ^' acquire vacuum 
expectation values, SU{2) x U{1) and Gp are broken and the neutrinos acquire masses given by the mass matrix 
M,j oc< ^ >< ^ > as well as the charged leptons. (Henceforth, for a Higgs doublet ^ we use < ^ > to denote the 
vacuum expectation of the lower electrically neutral component of ^.) 

Let Ml, be diagonalized by U'^M^Uv = so that the 3 neutrino fields that appear in are related to the neutrino 
fields v"^ with definite masses by = Ui.v'^ . Similarly, let the 3 charged left handed lepton fields I that appear in i/j^ 

iTjl'^^m I SO that the neutrino mixing 

[I^Oablb J 

matrix as defined in (1) is given hy V = UiU^. One difficulty in constructing a theory for V is that it arises from the 
"mismatch" between two rotations Ui and C/y. 
As in turns out, in our model building efforts, we often have to forbid the <^'s that appear in O4 from appearing in 

O5. This could easily be implemented by imposing a discrete symmetry under which ip e^^ip, 

iC ^ gixiC (where 

e*^ 7^ —1 is some appropriate phase factor), with all other fields unaffected. We will leave this implicit in what follows. 

Within the minimalist framework outlined here we offer some possible schemes. None of these could be said to be 
terribly compelling but at least we keep within the usual rules of the model building literature. The various schemes, 
depending on what representations of A4 we choose for the various fields tjj, , and ip, could be listed systematically. 

IV. Model A 



We first try the assignment V' ^ 3, 1, 1', and 1", and <^ ~ 3. The Lagrangian then contains the terms 

/ll^f ((/?|Vl + ¥'2V'2 + VsV's) + /l2^^(w(/5|V'l + 'y^'2^2 + W^t/JgV's) + h^l^ {t^"^ ^{i^l + <P2V'2 + ^(/sJV'a) (12) 

//ii \ /I 1 1 \ /v^I \ /v>i\ 
= ( /f ) h2 CO 1 p\ \\^^\ (13) 

V I \uj^ I u: I \ q pIJ J 

It is natural for the 3 < pa >= Wq's to be equal since A4 requires that the coefficients of pl^Pa and of (pl^pa)^ in the 
potential be independent of a = 1, 2, 3. (See the Appendix for a more detailed analysis.) If so, then upon spontaneous 
gauge symmetry breaking we obtain 
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{I? l§ l§)\ \ i CO 1 iv^ \ \ \ (14) 

\0 rur J \uj'^ 1 CO 

with rUe = hw and so on. It is useful to define the "magic" matrix [28] 

/I 11 \ 

A = \ CO 1 u;2 (15) 




Then = or / = y/SA-H"" so that U} = {y/3A-^)'' = -j=A 



The crucial observation at this point is that the sum of the first and third columns in A gives — 1 | and that 

V-1 

the difference of the first and third columns in A gives y/3i 11 1 , which up to some overall factors are precisely 
the first and third column respectively in the desired V in (5). In other words, if 



1 -1 

r- I 
v2 \ 1 1 



U, = ^\ V2 1 , (16) 



then U^Uiy = with V the desired mixing matrix in (5) and the diagonal phase matrix $ with the diagonal elements 
— 1, 1, and —i. Thus, if we could obtain [/^ we would achieve our goal of deriving V. 

We recognize that Ui, is just a rotation through 45° in the (1 — 3) plane. Recalling that Ui, is determined by 
requiring M^Uv = be diagonal we see that if we could obtain an of the form 



= 7 (17) 




(note that the 2x2 matrix in the (1 — 3) sector has equal diagonal elements) then we are done. Our discussion here 
overlaps with that given recently by Babu and He [22]; however, their discussion is given in the context of a much 
more elaborate scheme involving supersymmetry. 

Referring to (11) wc see that by imposing a discrete symmetry K2 under which -02 ^ ~'02, '^2 ^ ^^2, with all 
other fields unaffected, or equivalently a discrete symmetry i^is under which —^1,^3 —"03; 'Pi ~^ —fii 

ip3 — > —ip3, with all other fields unaffected, we can obtain the texture zeroes in (17), but unfortunately this does not 
imply that ( M^)ii = ( M^)^^. Furthermore, K13 is just the element r2 of A4 and so it does not commute with A4. 
Note that upon the ip^s acquiring equal vacumii expectation values, A^ is broken down to a Z3 generated by {/, c, a} 
and unfortunately r2 does not belong in Z3. Perhaps there is a more attractive scheme in which a reflection symmetry 
like K2 could emerge effectively. 

In another attempt to obtain an of the form in (17) wc introduce Higgs doublets x a-nd ^ transforming as I and 
3 respectively. We then have three types of O5 operators, namely (xi^)^, (xV')(?V')) and ■ As mentioned earlier, 
we impose a discrete symmetry to forbid (p from participating in O5. As discussed in the Appendix, we could naturally 
suppose that the vacuum expectation value of ^ points in the 2-direction, that is, < ^2 >7^ with < ^1 >=< ^3 >= 0. 
Let us now list how the different O5 operators contribute to M^, upon x and ^2 acquiring a vacuum expectation 
value. The operator (x'f/')^ contributes a term proportional to the identity matrix. Next, {x'^){^''P)^ which is formed 
by 3 X 3 X 3, consists of two terms, corresponding to the two ways of obtaining a 3 upon multiplying 3x3. One 
term has the form x{'^i^2ip3 + ^2C3V'i + V'3CiV'2), with the other term having an analogous form. Thus, the operator 
(xV')(^V') contributes the term denoted by /3 in (17). Finally, the operator (^V')^ actually denotes schematically 4 
different operators since it is formed by (3 x 3) x (3 x 3) and this contains 1 x 1, 1' x 1", 3 x 3, 3 x 3, and 3x3, 
corresponding respectively to the operators 

(6V'l+6V'2+6V'3)', (18) 

{^iipi + a;6V2 + w^^3V'3)(6V'i + <^^6V'2 + w^sV's), (19) 
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(20) 



(e3V'2, e2V'i)-(e3V'2, Civ's, ^v-o = (6V'2)' + (Civ's)' + (^V'O' (21) 

and 

(6V'3, 6^1, Cl^2).(6V'3, 6^1, 6^2) = (6^2)' + (6^3)' + (6^1)'. (22) 

(This is essentially the same as the analysis of an A4 invariant Higgs potential given in the Appendix.) Upon ^2 
acquiring a vacuum expectation value, we obtain respectively tjj2tp2, V'2V'2, 0, ViV'i and VsV's- Unfortunately, the 
effective coupling constants in front of the operator in (21) and (22) are in general not equal to each other and thus 
we obtain an Mi, of the form 

/ a-e 0/3 \ 
= 7 (23) 
\(3 a + e J 

rather than the M,y in (17). To set e to wo would have impose a discrete interchange symmetry P13 which interchanges 
the indices 1 and 3 but unfortunately, just as before for K13, P13 does not commute with A4. 

At this point, we could only suppose that s is small compared to /3, in which case Ui, is perturbed from the desired 
U„ in (16) to 

U. = ^\ V2 10 (24) 

V2 Vl 1 J \^ I J 

The resulting deviation from the V in (5) may be interesting phenomenologically. In particular, Ve3 — is no 

longer identically 0. In [11] it was advocated that experimental data be parametrized as a deviation from V in (5) as 

discussed in section III there. 

In this scheme, the neutrino masses come out to be a — y'/?^ + ,7, and a + \J (S^ ^ and thus both the normal 
hierarchy and the inverse hierarchy could be accomodated by suitable tuning, but there is no true understanding of 
neutrino masses as remarked earlier. 

V. Model B 

Following Ma [27], we take V" 3, Z*^ ~ 3, and ip 1,1', and 1". In other words, we have 3 Higgs doublets each 
transforming as a singlet under A^. The Lagrangian then contains the terms 

/ll</?I (Zf Vl + l2^2 + V-S^fV-s) + /l24(^f V'l + '^'??V^2 + t^V'S^fV's) + /i3<^k«f V-l + t^^FV'2 + W^i/.g^fi/.g) (25) 

Upon the (^'s acquiring vacuum expectation values v we obtain a diagonal charged lepton mass matrix, with the charged 
lepton masses given by the absolute values of h\V\-\- h^vi + /isiis, h\V\-\- uP'hiVi + ivhsvs, and hivi+ u}h2V2 + iv'^hsvs. 
All that matters here for our purposes is that we have enough freedom to match the observed masses me, m^, and 
rUr- The salient point here is that Ui — I, so that we only have to worry about getting the desired U,y. 

As is obvious and as was discussed in [11] and in [5], in a basis in which the charged lepton mass matrix is already 
diagonal, the neutrino mass matrix M^, is of course determined in terms of the three neutrino masses and the neutrino 
mixing matrix V. Call the three column vectors in the mixing matrix Uj. Then M^, is given by 

3 

= ^mm(tTi)^. (26) 

i=l 

In particular, if we believe in the V in (5) we have 

/lll\ _/oo 0\ 

(27) 

1 ly ^ \o -1 1 J 

1 1 i\ 

With A4 it is natural to obtain the matrix Mn = | 1 1 1 and the identity matrix. In particular, if we introduce 

1 1 ij 

a Higgs doublets C transforming as 3 under A4, and arrange the Higgs potential such that the 3 vacuum expectation 
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values < ^1 >=< ^2 >=< Cs > s-rc equal. We then see from the list of operators of the form {(.'p){£.f) given in (18-21) 
at the end of the last section that we obtain for an arbitrary linear combination of M]j and the identity matrix, 
which is not what we want. 

In [5], in discussing the neutrino mass matrix, we proposed a basis of 3 matrices other than those that appear in 
(27). First, the 3 column vectors in V are the eigenvectors of the matrix 



with eigenvalues mi = m2 = 2a, and ms = —4a. (The parameter a merely sets the overall scale.) Thus, with Mq 

as the mass matrix Arn^i = and this pattern reproduces the data | Ato2iI/I^™32I "C 1 to first approximation. 
Because of the degeneracy in the eigenvalue spectrum, V is not uniquely determined. To determine V, and at the 
same time to split the degeneracy between mi and m2, we perturb Mq to M = Mq + eaMo- The matrix Md is 
evidently a projection matrix that projects the first and third columns in V to zero. Thus, the eigenvalues are given 
by mi = 2a, 7712 = 2a(l + 3e/2), and = —4a, where to the lowest order e = H^m^i/ Am%2 and a? = Am32/12. 
Finally, to break the relation \m^\ = 2\mi\ ~ 2|m2| we can always add to M a term proportional to the identity 
matrix. But it seems difficult to get the matrix in (28) using A4 alone. 

VII. Other possibilities and conclusion 

Given that A4 has only 4 distinct representations we could of course systematically go through all possibilities. 
Thus, next we could take ~ 3, ^ 3, and ip ~ 1,1', and l". The charged lepton mass term would have a form 
analogous to that given in (12). But clearly, if we now assume the < >'s to be equal, we once again get the matrix 
A but now acting on 1*-^ instead of on ip. Note that if we assign ■i/'2 ~ 1 and ipi, ■03 to 1' and 1" respectively and 
introduce a Higgs doublets x transforming as 1, we get via the operator O5 a neutrino mass matrix M^, of the form 
in (17) but with a = 0. 

Another possibility is to assign V', (fi, and all to the 3 in which case the charged lepton mass matrix is generated 

by two terms, h{Lp\l2ip3 + (fill^'^Ji + (P3litp2) and h'{ip[l^'^2 + (P2^f V's + 'P3^2'V'i)- If we assume the < ^pa >'s to be 
equal, then the three charged lepton masses are given in terms of only two parameters. 

In conclusion, we have discusssed various schemes to obtain a particularly attractive neutrino mixing matrix that 
closely approximates the data. Instead of detailed models, wc use a low energy effective field theory approach, allowing 
only Higgs doublets to survive down to the electroweak scale. We have also explicitly made the restrictive assumption 
that A4, survives down to the SU{2) x U{1) breaking scale. Of course, if Higgs triplets could also be used, as for 
example in [16], or if A^ is broken at higher scale (for example by the coupling of the s^'s to the singlet scalar field h 
in the model in [26]), then many more possibilities open up and one could go beyond the discussion given here. We 
have been intentionally restrictive here. Ultimately, of course, any disucssion of neutrino mixing shuld be given in a 
grand unified framework (for recent attempts, see for example [29] [30] in which neutrino masses, as well as quark 
masses and mixing, are also "explained." We do not attempt this more ambitious program in this paper. 



We need to study the Higgs potential for several SU{2) x U{1) Higgs doublet i^'s which transform according to 
various representations under A4. For the sake of simplicity, here we restrict ourselves to the Higgs potential for a 
single SU{2) x C/(l) Higgs doublet ip which transform like a 3 under A4. Hopefully, the conclusions reached with this 
restricted analysis continue to hold when the couplings between different Higgs doublets are small. The multiplication 
3x3 = 1 + I'+l" + 3-1-3 tells us that there is only one quadratic invariant s = (flfi + ^If^ + 'filfs- 

Since ( 3 x 3) x (3 x 3) contains 1 four times, corresponding to 1 x 1, l' x l", 3 x 3, 3 x 3, and 3 x 3, we should have 
5 quartic invariants. The obvious quartic invariant is g = = ( ip\<fi + <fi2V2 + 'fit'fis)^ ■ Corresponding to V x V\ we 
have {(p\(pi +tJ(p\(p2 + ^^'^fi^ajifi'Pi +tj->^f2'P2 + '^'P^'Ps) giving rise to q and the quartic invariant q' — (p\Lpiipl(p2 + 
¥'2'^2¥'3<^3 + </'3<^3</'iVi- Next, corresponding to 3 X 3 and 3 X 3 we have g" = {(p\(p2, 'P2'P3, 'P3'Pi)-{'P2fi: 'P3'P2, 'Pi'fis) = 

bi</'2p + |<^2V'3l^ + and q'" = {ip\ip2,(plip3,V^3Vi)-iVi^2,V2^3,V3V>i) = {^1^2)^ + (Avs)'^ + ivWi)'^- The 

5th invariant is the complex conjugate of q'" . 

Thus, the most general Higgs potential is given hy V = —fj^s + Ag + A'g' + A"g" + \ (\"'q"' + h.c.) Assuming that 
the 3 (/7's all point in the same direction within SU{2), then we have V = —fj?'{v\ + + vf) + \{v\ + 1;* + vf) + 
\{v\v2 + vfi^f + vlvD where A = 2A + A' + A" + A'". For the sake of simplicity, we will take A'" and the various t;'s 
to be real, since our focus here is not on CP violation. 




(28) 



Appendix 
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It is then straightforward though tedious to calculate the value of V and the eigenvalues of the second derivative 
matrix evaluated at the three mimina of interest: E : {vi = V2 = V3 = v}, U : {vi = v,V2 = V3 = 0}, and 

P : {vi = V2 — v,V3 = 0}. We find 

C/ : = = -f„i^= [V, ^^^] and 

p . 1,2 _ y, iil, Q _ M 2 2^^(A-2A) 2/.^(A-2A) i 

~2A+A''^I^~ ^t^I' " - l*/^ ' 2A+A ' 2A+A J 

We note that by choosing A < and sufficiently close to —A or by not doing this we could set V\e much lower 
than V\u or vice versa. On the other hand, for P to be a minimum, we need A — 2A > 0, which would make V\u 
lower than V\p. It appears that in this simple one Higgs doublet case, P is never the true minimum. Of course, in all 
the models we discussed, we have to introduce more than one Higgs doublets and so presumably almost anything is 
possible by coupling the various doublets together. 
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